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Abstract. We study the phase diagram of two different Hamiltonians with 
compctiting local, nearest-neighbour, and mean-field couplings. The first example 
corresponds to the HMF Hamiltonian with an additional short-range interaction. 
The second example is a reduced Hamiltonian for dipolar layered spin structures, 
with a new feature with respect to the first example, the presence of anisotropics. 

The two examples are solved in both the canonical and the microcanonical 
ensemble using a combination of the min-max method with the transfer operator 
method. The phase diagrams present typical features of systems with long-range 
interactions: ensemble inequivalence, negative specific heat and temperature 
jumps. 

Moreover, in a given range of parameters, we report the signature of phase 
reentrance. This can also be interpreted as the presence of azeotropy with the 
creation of two first order phase transitions with ensemble inequivalence, as one 
parameter is varied continuously. 
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1. Introduction 

Interactions among elementary constituents of matter can be classified according 
to the features of the two-body potential. If this latter decays at large distance 
with an exponent that is bigger than space dimension, one speaks of short-range 
interactions: otherwise interactions are long-range [TJ [21 [31 H] • The difference between 
these two types of interactions has consequences at both the thermodynamic and 
dynamical level. In particular, systems with long-range interactions present a series 
of peculiar properties which are mainly a consequence of the lack of additivity of 
energy. Perhaps, the most impressive one is ensemble inequivalence [5], which entails 
the presence of negative specific heat in the microcanonical ensemble. Recently, 
there has been a growing interest in systems with long-range interactions. Besides 
realizing that they appear in many branches of physics (self-gravitating systems, two- 
dimensional hydrodynamics, dipolar interactions, unscreeened plasmas, etc.), the main 
breakthrough has been the discovery that simple toy models of the mean-field type 
reproduce many of the fundamental properties of systems with long-range interactions. 
These models can be solved exactly in both the canonical and the microcanonical 
ensemble, therefore ensemble inequivalence can be checked explicitly. Moreover, the 
dynamics of Hamiltonian systems with long-range interactions is well described by 
the Vlasov equation [B] , which allows one to interpret the presence of quasi-stationary 
states as stable stationary states of the Vlasov equation. 

An interesting but largely unexplored aspect of the physics of systems with long- 
range interactions is the combined effect of terms of short and long-range type. 
The Ising chain with mean-field (Curie- Weiss) interaction and nearest-neighbour 
interactions has been solved in the canonical ensemble by Nagle [7] and Kardar [5] (see 
also Ref. [9]). The complete microcanonical solution of this model has been recently 
obtained [TU] and it has been realized that the addition of a short-range term to the 
Curie- Weiss one is enough to generate ensemble inequivalence. The emphasis of this 
latter paper was rather on ergodicity breaking, a property which is induced by the lack 
of convexity of the accessible region of macroscopic thermodynamic parameters (here 
energy and magnetization) . It was later realized that ergodicity breaking can be found 
generically in systems with long-range interactions [Tj . 

After the introduction of an XY model with mean-field interactions [T2], the so- 
called HMF model, a model defined on a one-dimensional lattice and with additional 
nearest-neighbour interactions has been considered |13) . A preliminary study of the 
phase diagram and of the robustness of quasi-stationary states has been performed. As 
for the previous Ising-type model, ensemble inequivalence is induced by the addition of 
the short-range term, since inequivalence is absent in the pure mean-field XY model. 
In addition to results already published elsewhere [T31 H] , we present here an elegant 
method to derive the second order phase transition line using a perturbation theory 
inspired by quantum mechanics. 

There is a growing interest in finding physical systems in which it could be possible 
to observe experimentally the effects found in systems with long-range interactions. 
Examples are the free electron laser and the collective atomic recoil laser discussed in 
Ref. [14] in this issue, which should allow one to observe quasi-stationary states. 

For a certain class of layered magnets (e.g. (C 'H3N Hj,)2CuCIa) , it has been 
remarked [15 that, in specific temperature ranges and for given shapes of the samples, 
it is possible to reduce the microscopic Hamiltonian of the system to that of an effective 
Hamiltonian of classical rotators on a one-dimensional lattice with local, nearest- 
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neighbour and mean-field interactions. Indeed, this Hamiltonian turns out to be very 
similar to that introduced in Ref . [T3] . This reduction is generic for systems dominated 
by dipolar forces [161 117) . for which the energy is marginally superextensive. Indeed, 
it turns out that an extensive energy can always be obtained, but it contains a shape 
dependent "demagnetizing" term proportional to the square of the magnetization. In 
our case, this latter term has the form of a mean-field Curie- Weiss type interaction |18) . 

In this paper, we will determine the phase diagram of this effective Hamiltonian 
for different values of its parameters. We will also take the opportunity to explain in 
detail the method of solution in both the canonical and microcanonical ensemble. It 
consists of an interesting combination of the min-max method introduced in Ref. |19j 
with the transfer integral method [20l [2ll [22]. As for the transfer integral, we will 
apply the Bode discretization method [23], which allows us to obtain a much more 
precise integration of the transfer operator with respect to previous studies [13] . 

Moreover, a phenomenon of phase reentrance [241125) is observed for the model we 
discuss in this paper. When lowering the temperature, one observes a disorder/order 
transition followed by an order/disorder one, in such a way that the low temperature 
phase has zero magnetization. 

The paper is organized as follows. In Section [2] we introduce and show explicitly 
the solution of the XY model with short and long-range interactions [T3]- In Section[3j 
we study the phase diagram of the reduced Hamiltonian which describes the layered 
magnetic system and discuss the phenomenon of phase reentrance. In Section [4j we 
draw some conclusions. 

2. HMF model with additional short-range interactions 

Let us consider the Hamiltonian [IJ [T3] 

N 2 N N 

H = E f - K J2 cos - 6 ^ + 2N [1 " cos( ^ ~ 6j)] (1) 

i—1 i—1 i.j — 1 

where {6i,Pi) are canonical conjugate variables, with 9i representing an angle on the 
unit circle. This Hamiltonian models a ID lattice whose sites are occupied by rotors: 
the lattice has periodic boundary conditions 9i+N — &i- Two kinds of interaction are 
present: one with coupling K, which is nearest-neighbour, and one with ferromagnetic 
coupling J, which is of mean-field type. We will allow the coupling K to take both 
positive and negative values. When K — 0, one recovers the Hamiltonian Mean 
Field (HMF) model, whose phase diagram in both the canonical and microcanonical 
ensemble has been widely studied. The model displays a second order phase transition 
at the critical temperature T c — J/2 and the corresponding critical energy u c = 3 J/4. 
Conversely, when J = 0, the model reduces to the well known chain of coupled rotators, 
which may have interesting dynamical effects [25] but does not have phase transitions. 
When combined, the two interaction terms produce a very rich phase diagram with 
first and second order phase transitions and a tricritical point. For convenience, since 
we will always use ferromagnetic mean-field interactions all along this paper, such as 
to generate order in ID, we will set J = 1. It has been shown [T3] that ensembles are 
not equivalent for this model and that the canonical phase diagram differs from the 
microcanonical one. All the interesting features of models with long-range interactions 
are present: negative specific heat, temperature jumps, breaking of ergodicity £Q. 
Here, we will not review these results, but we will introduce the method which allows 
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to solve the model, because it will be used also for the variant of model ([I]) that 
we present in this paper. Our theoretical approach combines three methods: one is 
specific of short-range interactions, the transfer integral method [201 US H2]; the second 
is standard for mean-field models, the saddle point method |29j : the third, the min-max 
method [H] , has been more recently introduced and allows one to get microcanonical 
entropy from the canonical free energy. Other solution methods introduced for long- 
range interactions, like the large-deviation method [30[|3T], cannot be simply extended 
to this class of model because of the short-range component of the Hamiltonian, which 
does not allow one to write the Hamiltonian in terms of "global" variables pQ. 

The phase transitions we will discuss in this paper are all characterized by the 
XY magnetization 
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which is a vector whose modulus m is finite in the broken symmetry phase. 

In order to derive the equilibrium properties of model , one first computes the 
partition function in the canonical ensemble as follows 



N 

Z = J n( d ^)exp(-/3H) 



(3) 



2tt\ 2 / j3N 
exp 



A? 



N 



exp 



i=l 



^sinfli 



2N 



N 



N 



2^ cos Qi 



\i=l 



W\- (4) 



Since two quadratic terms are present in the Hamiltonian, one has to introduce two 
Hubbard-Stratonovich transformations 
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with auxiliary fields z\ and z 2 , obtaining 
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where Zi — j3zi. It can be easily shown that, due to invariance under rotation in 0j, 
only the integral in the modulus z — \J z\ + z\ appears leading to 

2tt\ 2 / /3N\ fiN " 



( 0N\ fiN f n T-r 

exp I — I — / 2-kzAz Y | aOi exp 




i=i 

N N 

+f3zJ2cos8 l + /3kY^ cos (6 i+1 - 0;) ) . (8) 
i=i »=i 
The integral over the angles has the following form 

JV / N \ 

JJd^exp [Y,9{0i,0 i+ i)) , (9) 



where 



g (e, 9') = i/3z(cos 6 + cos 6') + (3K cos{0 - 9') . (10) 

Using the integral operator representation 

(Ti/>)(9)= [ d6'exp[g(9,e')]TP(6') (11) 

and taking advantage of the periodic boundary conditions 6\ = 9 N+ i, integral Q 
reduces to [2"2"] 

oo 

E^f (12) 

3=1 



where \j is the discrete set of eigenvalues of the transfer operator (11). When N is 
large, the maximal eigenvalue A dominates in the sum (corrections being exponentially 
small in N). Since the transfer operator is symmetric, all eigenvalues are real. 
Moreover, we expect on a physical basis that the maximal eigenvalue be positive, 
because a negative eigenvalue would generate an imaginary part of the free energy 
(see below) and would produce an oscillation in N spoiling the convergence to the 
thermodynamic limit. Expressed in terms of the maximal eigenvalue A of the transfer 
operator, the partition function reads 

2n \^-„f @ N \ „ /"..j.. _( N/3 _ 2 



Z= I j exp^-^-J/3iVyzdzexp^-^z 2 + 7VlnA(/3z,^)J .(13) 

Now, one can compute the partition function using the saddle-point method to perform 
the integral in z in the large N limit. The z factor in front of the integral gives a 
contribution which is negligible in this limit. One gets the free energy /(/?) as 

-Pf(/3)= -<p{0)= lim ilnZ (14) 

N^-oo iv 



sup 



^H±^l +lnX (f3z,Kf3)+ 1 -\n 2 ^ 



(15) 



= sup [-^(,3, z)] (16) 

z 

where the rescaled free energy <p = /3f has been introduced for convenience, because it 
appears naturally in the applicatiom of the min-max procedure Q]5]. In the second line, 
we also define the rescaled free energy ip as a function of both the inverse temperature 
and the Hubbard-Stratonovich auxiliary variable z. 
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One obtains the microcanonical entropy using the min-max method |19j 

s(u) = supinf[/3u + tp(/3, z)] (17) 

z ft 



sup inf 

/ ft 



0u - /?^^ + In X(6z, K(3) + 



(18) 

sup [s(m, z)] , (19) 



where s is defined in analogy with (p. Solving the variational problem in formula (18) 
also gives the value of the microcanonical temperature as a function of energy u. 

The second order phase transition line in both the canonical and the 
microcanonical ensemble can be derived analytically using a perturbation theory for 
the transfer operator inspired by quantum mechanics. This is possible because the 
magnetization m varies continuously at the phase transition, where it vanishes. Hence 
all quantities can be developed in Taylor series around z = 0, e.g. the free energy, 
which is even in z, reads 

-<fi(0,z) = a + bz 2 + cz A + o(z 6 ) . (20) 

The change of sign of the z 2 term b (with c remaining positive) will indicate where 
the second order phase transition is located in the (K, T) phase plane. The canonical 
tricritical point is obtained by requiring that also c vanishes. Therefore, in order to 
determine the second order phase transition line in the (if, T) plane, it is enough to 



develop the maximal eigenvalue of the transfer operator, appearing in (14), to second 
order in z. The series development of the operator is 

(Ttf)(0) = £(2WO(0y , (21) 
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Since we have to perform the calculation up to second order in z, we need to consider 
the perturbation theory for the transfer operator up to second order. We will 
concentrate on the perturbative corrections to the maximal eigenvalue, considering 
as the unperturbed problem the one with respect to the operator To, which can be 
solved exactly (see below). Let \r) be the eigenvector of To with eigenvalue A° 

T \r) = X° r \r) , (23) 

the maximal eigenvalue and eigenvector corresponding to the label r = 0. Non 
degenerate perturbation theory for Hermitian operators tells us that 



A = A" + zTOlO) + z 2 | (0|T 2 |0) + J2 'a™' I + ° (z4) • (24) 

It is well known that, introducing I r is the modified Bcsscl functions of order r, the 

eigenvalues of To are 2irI r (0K) with eigenvector exp(ir9) / V2tt (plane waves). It is 
then straightforward to compute the different terms of the perturbation series 

(0|li|0) = (25) 

(r\Tx\0) = ^(S rA + £ r ,_i) + hifiK)) (26) 

(0|T 2 |0) = ^ 7r (/ (^)+/i( i 9if)). (27) 
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Finally, one obtains the following formula for the maximal eigenvalue up to order z 2 

\(Rv RK\ n —T (RK\ i (P*) 2 *MPK) ( I (f3K) + h(f3K) \ 
\{Pz,PK) = 2irI (PK) + { hm) _ hm) ) 

= A + Bz 2 + ...., (29) 
from which we can compute the coefficients a and b of the free energy expansion in z 

a= -^ + \^J+^ A (30) 

By requiring 6 = 0, one gets the second order phase transition line 
h(pK) 2-/3 



h{PK) 2 + /?• (32) 
Indeed, it is an implicit equation linking j3 and K. To compute the following order 
in the free energy, z 4 , and the corresponding coefficient c, it would be necessary 
to perform calculations to second order in perturbation theory for the transfer 
operator. This could in principle be easily done, but the calculations are lengthy. 
In order to locate the canonical tricritical point, we have preferred to rely upon a 
numerical method. Since the coefficients a and b are determined analytically, one 
can numerically check the difference <p(/3, z) — a — bz 2 and observe when, for fixed K 
and varying T, this difference turns from positive to negative. In this way we get 
K CT p - -0.171, T CT p - 0.267. 

The first order phase transition line is determined numerically by finding, for each 
value of K, the temperature T t — l//3 t at which the value of the free energy at the 
two minima in z — and z^O coincide, i.e. <p(0t, 0) = z), where z is the jump 
in magnetization at the phase transition. 

One point of this line, the one at T = 0, can be computed analytically. Indeed, 
at T = the free energy coincides with the energy and, therefore, one can locate the 
transition value of K by requiring that the energy in the two phases coincides. At 
T = and negative K, one has a transition from a ferromagnetic phase, where the 
mean-field term of the Hamiltonian dominates over the short-range antiferromagnetic 
one, and an antiferromagnetic phase governed by the short-range contribution. One 
has just to choose typical microstates corresponding to the two phases and compute 
their energy. The kinetic energy can be taken at zero, since we are at zero temperature. 
In the ferromagnetic phase, the magnetization vector is invariant under rotation; wc 
can therefore assume that all the angles are at 0, = 0, giving the energy 

[/> = - KN - ^(Nf - ^(0) 2 = -n(k+^J, (33) 

disregarding constant terms. On the contrary, in the antiferromagnetic phase the 
angles are alternatively disposed at and n along the lattice, leading to the energy 

U< =0 + KN-^(0) 2 -^(0) 2 = +NK . (34) 

Imposing that [/> = U < , one gets the transition value of K: K trans = —1/4. 

The phase diagrams in both the canonical and microcanonical ensembles can be 
directly derived from the properties of the functions — <p(/3, z) and s(u, z). 

In the canonical ensemble, the free energy is determined by a single maximization 
procedure over z. After fixing the value of K, one looks for the value of (3 at which 
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the maximum changes from z = to z ^ 0: this can happen in a discontinuous way 
for a first order phase transition or continuously for a second order phase transition. 
Indeed, when the transition is first order the local maximum of — (p{P, z) at z = 
remains a local maximum. For second order phase transitions the local maximum at 
z = becomes a minimum. The tricritical point is found by determining when one 
changes from the first type of behavior to the other: in practice one observes that, 
coming from the first order phase transition side, a local minimum in z ^ merges 
with the local maximum in z = and converts it into a local minimum. 

Numerically, the crucial point is the determination of the maximal eigenvalue 
A(/3z, A/3) of the transfer operator. The integral in Eq. (11) is discretized using 
the Bode method 23 , which turns out to be much more efficient than the more 
standard trapezoidal discretization, allowing us the calculation of the eigenvalue with 
high precision using about 100 discretization points (to obtain the same precision 
with the trapezoidal method one should have worked with about 1000 points). The 
calculation of the maximal eigenvalue of the, now discrete, transfer operator is done 
using a routine of the Lapack library [3"2"| . 

In the microcanonical ensemble, the procedure is more complex because it requires 
the solution of a double extremalization, see Eq. ( 18 ). Moreover, one has to choose the 
appropriate range of energies u around the phase transition. Hence, one first fixes the 
value of u and performs a minimization with respect to (3, and then finds a function 
of z. The maximum of this latter function determines the equilibrium magnetization. 
The search for the transition lines of first and second order is in this respect similar 
to what has been already described for the canonical ensemble. The only important 
variant is that for microcanonical first order phase transitions, there is a temperature 
jump [5 at the transition energy. The two temperatures at the jump are found from 
the two different values of the minimum in j3 at z — and z ^ when solving the 
variational problem in /3 in Eq. (18). Practically, the microcanonical solution requires 
an order of magnitude more of computational effort, which is due to the need of 
diagonalizing the transfer integral more often. 

The phase diagram in the (A, T) plane is represented in Fig. [I] For positive 
values of A there is no competition between the short-range and the XY mean-field 
term in the Hamiltonian, and therefore one expects, and finds, a second order phase 
transition in both the canonical and the microcanonical ensemble (thick full line). 
For negative values of A there is instead competition and one can predict that for a 
sufficiently negative value of A, the phase transition becomes first order. Indeed, we 
have already given a theoretical argument, valid at T = 0, that suggests a first-order 
phase transition at K tra ns — —1/4. Hence, we can anticipate the existence of a first 
order phase transition line in the (A, T) plane and the presence of a tricritical point. 
What we find is analogous to what has been also derived for the Kardar-Nagle [8] 
Ising spin chain in one dimension |10j . It is in the range of A where the transition 
is first-order that ensemble inequivalence takes place. As it happens also for other 
models [T], the microcanonical (MTP) and canonical (CTP) tricritical points do not 
coincide. In the canonical ensemble, the transition line changes to first order (dashed 
line) at the canonical tricritical point (CTP), located at (A" = —0.171, T = 0.267) 
and ends at T = when A = K trans = —1/4. In the microcanonical ensemble, the 
transition is also second order on the right side of the diagram but remains second order 
also beyond the CTP (thin full line) in a region where the canonical ensemble predicts 
a first order phase transition. The second order line ends at the microcanonical 
tricritical point (MTP), located at (A = —0.182, T = 0.234). Beyond the MTP 
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the micro-canonical transition splits into two lines (dotted lines), corresponding to the 
two temperatures present at the transition energy (this is the so-called temperature 
jump); these lines finally merge at the zero temperature transition point. The two lines 
define a region which is not accessible in the microcanonical ensemble. Indeed, the 
microcanonical first-order transition corresponds to a single line in the (K, u) plane; 
it splits in the (K, T) plane because different values of T are found approaching the 
transition energy from above or from below. 
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Figure 1. Phase diagram of model Q in the (K, T) plane for both the canonical 
and the microcanonical ensemble. The second order phase transition in the 
canonical ensemble is represented by the thick full line. The transition changes 
order at the canonical tricritical point (CTP), shown by a filled triangle, becoming 
first order (dashed line). In the microcanonical ensemble, the transition remains 
second order (thin full line) also beyond the canonical tricritical point, reaching 
the microcanonical tricritical point (MTP). The first order phase transition in the 
microcanonical ensemble is represented by the two dotted lines, as detailed in the 
text. 

In summary, this simple spin chain system with short and long-range interactions 
is another example of system for which ensembles can be inequivalent. This means that 
experiments realized in isolated systems, described by the microcanonical ensemble, 
may give different results from similar experiments performed with well thermalized 
systems, for which the canonical ensemble is the appropriate one. It is interesting 
to emphasize that, although ensembles are equivalent for systems with short-range 
interactions, it is the introduction of the short-range term which has induced the 
inequivalence, since inequivalence is absent in the pure mean-field XY model. 



Models with short and long-range interactions: phase diagram and reentrant phase 10 



3. A reduced model for layered spin structures 

Spin systems with both Heisenberg short-range terms and dipolar interactions have 
been studied by Sievers and coworkers |15j . The spins are organized on N layers 
and the interlayer distance is sufficiently large that dipolar interaction dominates 
over Heisenberg exchange interaction. On the contrary the intralayer interaction is 
predominantly Heisenberg and induces a ferromagnetic coupling between neighbouring 
spins within a layer. In Ref. [18] it has been shown that, for rod-shaped samples, the 
Hamiltonian describing the system can be effectively reduced to the following one 

N 2 N B N 

ff= El^E cos -«i) + yE sin2 6 i 




(35) 

where 6j is an angle representing the orientation of a classical spin which describes 
the macroscopic magnetization in the j-th layer (the modulus of such spin is constant 
and can be set to unity). The constants K and B are determined in Ref. 18 from 
the microscopic parameters of the sample. 

Hamiltonian (35) looks very close to the one of Eq. ([!]). However, there are 



two substantial differences: i) an asymmetry in the mean field term, with the term 
containing the sines of the angles which is now positive, ii) the addition of a local 
potential proportional to B. Both of these terms break the rotational invariance of the 
Hamiltonian and make the behavior of the system richer. Physically, the asymmetry 
in the mean field term, which favours the orientation of the magnetization vector 
along the 9j — axis (the easy-axis along the rod) is originated by the presence of 
an antiferromagnetic "demagnetizing term" proportional to the magnetization square 
of the sample. The tendency to align the spins along 9j — is strengthened by the 
local potential term when B > 0, and is instead opposed by the local term when 
B < 0. Also the role of the nearest-neighbour term is important, as we have seen 
for Hamiltonian 0, since, when K < 0, it favours antiferromagnetic configurations, 
without however fixing the global orientation of the spins. 

After performing the Hubbard-Stratonovich transformation, the partition 
function reads 



i=l 




N 

exp ( - zf - Nflz'z + /3z i y ] cos 9 t +i/3z 2 sin( 



+f3K}^ cos (9 l+1 - 9i) -'— ^ sinO] , (36) 

8=1 t=l / 

and the g(9, 9') function entering in the transfer matrix becomes 

g{6,6') = -/3zi(cos0 + cosO') + ^f3z 2 (sm9 + sin 6') + f3Kcos(9 - 9') 

-^(sm 2 9 + sm 2 9'). (37) 
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The free-energy can be obtained, as before, by solving the following variational 
problem 

-pf(j3) = -¥>(/3)= lim ilnZ (38) 

N— foo iv 



sup 

21,22 



(39) 



2 /3 

sup [-<p(P,zi,zz)] , (40) 



2l,2 2 



where A is the maximal eigenvalue of the transfer operator. Analogously to Eq. ( |18[ ), 
one defines a microcanonical variational problem. 

There are two additional difficulties, when determining the free energy of this 
model. The first one is that, due to the lack of rotational symmetry, one has to look for 
extrema in two coordinates, z\ and zi. The second one is that, due to the positive sign 
in one of the mean-field terms, the function g(9, 9') is complex and, consequently, the 
transfer operator is no more Hermitian. This latter property entails that the maximal 
eigenvalue could have an imaginary part. However, this is not possible, because we 
have started our calculation from a partition function which is real, and in no way it 
can become complex in the end. We have then checked numerically that, indeed, the 
imaginary part of the maximal eigenvalue of the transfer matrix is by twelve orders 
of magnitude smaller than the real part, and that its size decreases further as the 
number of points in the discretization of the integral by the Bode method increases. 
We can thus safely assume that the maximal eigenvalue is real. There is moreover one 
further simplification; by the numerical search of the extremum for different values of 
the parameters, always for T > (we will come back to this point), we have found that 
the non trivial maximum of —(p(J3, z\, z%) always appears at z\ ^ and Z2 — 0. This 
means that the symmetry is always broken is such a way that the total magnetization 
is directed along the 9 = axis. This is of course perfectly understandable when 
B > 0. However, when B < the local term will tend to orient the magnetization 
at 9 = 7r/2 . What we find empirically then means that the mean-field term always 
dominates at T > and determines an orientation of the magnetization along the 
9 = axis. One can thus restrict the search of the extrema of — z\, z%) to the 



region where Z2 = 0, which also implies that the imaginary term in (37) vanishes, 
simplifying also the calculation of the maximal eigenvalue of the transfer operator, 
which is now necessarily real. 

The phase diagram for B > does not show any new features with respect to 
the one of Hamiltonian ([lj. This can be intuitively explained by the fact that the 
positivity of both the local term and the mean-field term with the sines of the angles 



in Hamiltonian (35) requires, by energy minimization, that magnetization aligns along 
the 9 = axis. This exactly reduces Hamiltonian ( 35 1 to ([T]) . This does not mean 



that one gets the transition temperatures which will be determined also by the effect 
of the entropy. Indeed, at fixed K, transition temperatures increase as B is increased, 
and the tricritical points move towards higher values of K. However, qualitatively, 
we do not expect significant changes of the phase diagram. Moreover, the transition 
at T = happens again at K trans = —1/4 for all values of B > 0, as can be easily 
proven by repeating the argument given in Section [2j In Fig.[2j we display the phase 



diagram of Hamiltonian (35) for B = 0.5 in both the canonical and the microcanonical 
ensemble. The transition lines and the tricritical points are drawn with the same line 
types and symbols as those of Fig. [I] In this positive B domain, the phase diagram is 
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essentially similar to what was obtained for model (FTl). 
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Eh 0.2 
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Figure 2. Phase diagram of Hamiltonian ( |35| for B = 0.5. The symbols and line 
types are the same as those of Fig. ^ The location of the canonical (CTP) and 
microcanonical tricritical point (MTP) are at (K = —0.1805, T = 0.2913) (filled 
triangle), (K = —0.1656, T = 0.3319) (filled circle), respectively. 



The case B < is instead more interesting. One observes numerically that, not 
only the transition temperatures, as expected for reason of symmetry with respect to 
the previous B > case, reduce (at fixed K) as B becomes more negative, but that 
the transition lines bend, with the lower part of the transition lines themselves which 
slides towards larger values of K. 

The transition at T — is also shifted towards a larger value of K . Can this be 
understood? The only way this can be justified is that the local term, proportional 
to B, plays a role. This can happen, at T = 0, only if the antiferromagnetic zero- 
magnetization state has its spins oriented along the 9 — ir/2 axis. Let us then repeat 
the zero-temperature argument of the previous section for the sake of clarity. In the 
magnetized phase, the mean-field terms dominate and the spins will be all oriented at 
6i = or it. The corresponding energy will be 



{/> = — KN 



1 

2N 



(AO 2 



1 

47V 



B 



(0) 



-N K 



(41) 



In the zero-magnetization phase instead, the local and nearest-neighbour terms in the 
Hamiltonian will dominate, and the microstates will be antiferromagnetic, with the 
spins alternatively oriented at 7r/2 and —tt/2. This gives an energy 



E7< = + KN 



V ) 2 + ^7(°) 2 + ?W- 



T ' 4N" ' 2 
Equating the two energies, U > — U < , one obtains K tra ns = —{B 
is in perfect agreement with the one determined numerically. 



l)/4. 



(42) 
This value 
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Figure 3. Phase diagram of Hamiltonian | |35| l for B = —0.5. The symbols and 
line types are the same as those of Fig. [I] The location of the canonical (CTP) 
and microcanonical tricritical point (MTP) are at (K = — 0.1623, T = 0.1605) 
(filled triangle), (K = -0.1626, T = 0.1574) (filled circle), respectively. 
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Figure 4. Zoom on figure [3] around the canonical and microcanonical tricritical 
points. 
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For B sufficiently negative, e.g. B = —0.5 as in Fig. |3j one observes the 
curious phenomenon of phase reentrance: a vertical line in the (K, T) plane crosses 
twice the transition line. By lowering the temperature at a fixed value of K, 
one observes first a disorder/order transition and then, by further lowering the 
temperature an order/disorder transition. Fig. [4] presents a zoom of this interesting 
region. Interestingly, the microcanonical first-order phase transition lines (dotted) 
seem to coalesce at the bending point before separating again below it. The 
system is disordered down to zero temperature. However, in our case, as we have 
commented above, the microstates close to zero temperature, although they have a 
zero magnetization, are antiferromagnetic states. 

Phase reentrance was first described by Griffiths and Wheeler [24] and has found 
later several experimental verifications for colloids and polymers. Recently, it has 
been also discovered in the context of out-of-equilibrium phase transitions of the HMF 
model [55]. 

A different interpretation of the phenomenon of phase transition line bending that 
we observe can be given in terms of the classification of phase transitions for long-range 
systems described in Ref. [26]. Let us refer for example to the situation represented 
in Fig. [4] When increasing K from negative values one passes from a situation where 
no phase transition is present (one is always in the m — phase for all temperatures) 
to one where two first-order phase transitions are present in both the canonical and 
microcanonical ensemble. This is a case of azeotropy with ensemble inequivalence 
discussed in Ref. [26]. To our knowledge, this kind of azeotropy has never been found. 
A different case of azeotropy at a second-order phase transition, i.e. with ensemble 
equivalence, was indeed found for certain two-dimensional flows in Ref. [27] . This 
type of phase transition pattern could be present in our model for smaller values of 
B, when presumably both the canonical and microcanonical tricritical point lie below 
the bending point of the transition line. Moreover, a case of second-order azeotropy 
with ensemble inequivalence could be found when the canonical tricritical point is 
above the bending point and the microcanonical one is below. All this remains to be 
ascertained. 

Among the curious effects related to phase reentrance is the fact that, even 
in mean-field models, one can have non classical exponents at second order phase 
transitions. This would happen in our model in the latter case we have discussed in 
connection with second-order azeotropy, in which the value of K is tuned in such a 
way that a vertical line in the (K, T) plane touches tangentially the transition line of 
second order. 

The phase reentrance phenomenon is, in our model, more and more pronounced as 
B becomes more negative. However, if B is so negative that the phase transition point 
at T = moves to positive K, we expect a change in the nature of the phase transition, 
because the microstate cannot be antiferromagnetic in this case. By varying the value 
of B, different patterns of phase transitions can be observed. Increasing temperature 
one can have a first order phase transition followed by a second order one, as shown 
in Fig. [5j Otherwise, one could first meet a first order phase transition and then a 
second order one or, finally, meet two second order phase transitions. 

4. Conclusions 

In this paper we have studied in full details the phase diagram of a Hamiltonian with 
local, nearest-neighbour and mean-field couplings which is obtained from a microscopic 
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Figure 5. Caloric curve and magnetization vs. energy at B = 0.5 and K = —0.15. 
Microcanonical temperature (full line) is plotted vs. energy (left axis) showing 
first a temperature jump, corresponding to a first order microcanonical phase 
transition, followed by second order phase transition detected by a change in 
slope. Magnetization values (dashed line) are shown on the right axis. 



model of dipolar layered magnets, e.g. {CH^NH^iCuCli [151 118] . 

Such layered magnets share the feature that the spins (here the Cu atoms) are 
organized on planes, and below a given temperature they order ferromagnetically 
on the plane. However, since the distance between layers is such that the Heisenberg 
exchange interaction becomes negligible with respect to dipolar interactions, the global 
ferromagnetic order is reached at much lower temperatures. For a sample of ellipsoidal 
form, the dipolar interaction gives rise to a term in the energy which is proportional 
to the square of the total magnetization, with a coupling constant which depends on 
the geometric features of the sample. This heuristically explains why in our effective 
Hamiltonian ( 35 1 , we have both a term describing the interaction between nearest- 
neighbour "classical" XY spins, which represents the magnetization in a layer, and 
mean-field terms which take the dipolar interaction between all the spins in the sample 
into account. 

Ising models with both short-range and long-range terms in the Hamiltonian 
have been independently considered HI [TU] and studied in both the canonical and 
microcanonical ensemble. They show all the typical features of systems with long- 
range interactions [H El El H], i.e. ensemble inequivalence, negative specific heat, 
temperature jumps, breaking of ergodicity and quasi-stationary states. Predictions in 
the microcanonical ensemble could be checked in experiments performed in conditions 
of thermal isolation or fast variation of external parameters (temperature quenches, 
external magnetic field sweeps) [33] . 

Besides the Ising models mentioned above, XY models on one-dimensional 
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lattices with both short and mean-field terms have been studied |13j . A preliminar 
reconstruction of the phase diagram and a study of the dynamics have been performed. 

With the aim of looking for experimental realizations of long-range effects, in 
Ref. [IS] an attempt has been made to fill the gap between the Hamiltonians introduced 
for magnetic layers dipolar systems and the XY Hamiltonians studied in Ref. [13 . This 
has led us to study Hamiltonians of the type of Eq. ( 35 1 . 

A new feature of this Hamiltonian, with respect to the previously studied XY 
Hamiltonian, is the presence of anisotropies. The first one is due to a different coupling 
(even in sign) of mean-field terms in cosines and sines of the angles of the XY rotor. 
The second is due to a local term, proportional to the coupling B, analogous to an 
external magnetic field. The value of this coupling constant B is calculable from the 
geometric and physico-chemical properties of the sample and turns out to be positive. 
However, in order to perform the analysis in full generality, we have considered both 
positive and negative values of B. 

The solution of Hamiltonian ( |35| is possible via a smart combination of the min- 
max method introduced in Ref. [ill] with the transfer operator method [SU1 |2"T1 |2"2"] . 
The diagonalization of the transfer operator has been performed using the Bode 
discretization 23J, which allows to work with substantially smaller matrices with 
respect to the more standard trapezoidal discretization. These methods allowed us to 
obtain a greater precision in the determination of the phase diagram with respect to 
previous studies of Hamiltonians of this kind. 

We have found that for B > the phase diagram is qualitatively similar to the 
one of the XY Hamiltonian ([!]). In the canonical ensemble a second order phase 
transition is found when the nearest-neighbour coupling K is positive. The phase 
transition becomes first order if K is sufficiently negative, inducing antiferromagnetic 
patches in the microscopic spin configurations (although antiferromagnetic global 
order is not possible in one dimension). The second order phase transition line 
is separated from the first order one by a tricritical point. In the microcanonical 
ensemble, the second order phase transition line coincides with the one of the canonical 
ensemble until the canonical tricritical point is reached. Beyond this point, the 
phase transition remains second order in the microcanonical ensemble until the line 
reaches a microcanonical tricritical point, which differs from the canonical one. At 
the microcanonical tricritical point the transition line splits in two, because of the 
presence of a first-order microcanonical transition. The two lines join again at zero 
temperature. These features are common with those found for the phase diagram 
of other systems with long-range interactions and are an impressive manifestation of 
ensemble inequivalence. 

Definitely more interesting is the B < case. The phase transition lines bend 
at lower temperatures as B becomes more negative until the lines take a boomerang 
shape. This is a signature of the phenomenon of phase reentrance, first discussed by 
Griffiths and Wheeler [24J. When lowering the temperature at fixed K the system 
first go through a disorder/order transition and then again through an order/disorder 
one and is, counter intuitively, disordered down to zero temperature. An alternative 
interpretation of this phenomenon makes reference to azeotropy \20>\ I27j . The two 
transitions can be of first or second order depending on the chosen value of B. In 
the microcanonical ensemble, also the splitted first-order phase transition lines bend. 
This latter effect determines the presence of a different caloric curve which shows a 
temperature jump followed by a more standard change in slope at the second order 
phase transition. 
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Much remains to be done in the study of this kind of models. One could mention 
that the solution of a more complicated Hamiltonian, which takes into account also the 
off-plane motion of the spins [TB], is possible using the same methods. The additional 
difficulty is that one has to diagonalize a transfer operator containing two angles and 
this makes the application of standard numerical methods quite hard. 

Fully absent in this work is the study of dynamical effects. Both breaking of 
ergodicity and quasi-stationary states should be present in Hamiltonian ( 35 ) and would 
be certainly interesting to investigate. 
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